Quantum Dynamics of a Bose Superfluid Vortex 
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We derive a fully quantum-mechanical equation of motion for a vortex in a 2-dimensional Bose 
superfluid in the temperature regime where the normal fluid density Pn{T) is small. The coupling 
between the vortex "zero mode" and the quasiparticles has no term linear in the quasiparticle 
variables - the lowest-order coupling is quadratic. We find that as a function of the dimensionless 
frequency fl — HQ/ItBT, the standard Hall-Vinen-Iordanskii equations are valid when S7 <C 1 (the 
"classical regime"), but elsewhere, the equations of motion become highly retarded, with significant 
experimental implications when $7 > 1. 
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Quantum vortices were first predicted in '*He superfluid 
by Onsager [1] and found experimentally a decade later 
[5]. Vortices, along with quasiparticles, constitute the 
two basic excitations in many condensed matter systems 
[5]; they may also have existed as topological defects in 
the early universe @]. 

Remarkably, the fundamental question of how vortices 
and quasiparticles interact, and how vortices move is very 
controversial, notably for superfluids The argu- 

ment is usually phrased in terms of forces acting on a 
vortex moving with respect to both the superfluid [hav- 
ing local velocity Vs(r) and superfluid density ps] and the 
normal fluid [with velocity f„(r) and density Then, 
the classical equation of motion for the semiclassical vor- 
tex coordinate Rv(t) (here taken to be a point in the 
plane - we discuss the 3D problem at the end) is usually 
written [8] as 

M,R, - fM - fqp - Fac{t) = (1) 

where Fadt) is some driving force. My is the vortex mass, 
/a/ = PsI^ ><■ (R-u — Vs) is the (uncontroversial) Magnus 
force for a vortex with circulation k = zh/m, and the 
quasiparticle force is 

fqp = Do{Vn - R,,) + D'^i X {Vn - R„) (2) 

where Do{T), D'^T) depend strongly on the temperature 
T. The classic discussion of lordanskii yields 

D'ST) - -np,,{T) (3) 

These "HVI equations," due to Hall, Vinen, and lor- 
danskii [H [9] , have been used to analyze thousands of ex- 
periments in superfluids and superconductors in the last 
60 years [TD]. However, there is no consensus on the value 
of either the mass My (estimates in the literature range 
from zero to inflnity (Tj) or of the coefficients Do{T), 
D'^{T). Indeed, Thouless et. al. ^ find D'^{T) = 
for all T, and scattering analyses give various results for 
D'^{T) [i ini [ITHI3]- The resolution of these questions 



has become an important unsolved problem in physics. 
We briefly discuss the experimental situation below. 

Previous analyses have been restricted to a local (in 
space and time) description, derived from semiclassical 
or perturbative calculations of quasiparticle scattering off 
a static vortex, with no vortex recoil [51 IHl ITT| - [T5] . This 
yields forces acting instantaneously on a quasiclassical 
vortex coordinate R^(i). There is no general agreement 
between these calculations (which are rendered difficult 
by the long-range vortex profile). Our tactic has been 
to formulate the problem fully quantum mechanically, in 
terms of an equation of motion for the vortex reduced 
density matrix p(R, R'; t) which is obtained by integrat- 
ing out all quasiparticle degrees of freedom. Here the 
vortex coordinate states |R), |R') are defined by the po- 
sition of the vortex node appearing in the many-body 
wave function. We then define a vortex "center of mass" 
coordinate R„ = ^(R-|-R'), and a "quantum fluctuation 
coordinate" ^ = R — R'. Equations of motion are then 
found for these two quantum variables, with the vortex 
recoil automatically incorporated. Remarkably, in ther- 
mal equilibrium we flnd that the results largely depend 
on one key parameter, the ratio SI = Ml/kBT, where f2 
is the characteristic frequency of the vortex motion, and 
fcsT is the thermal energy of the quasiparticles. When 
O <C 1 we are in a "classical regime," where we find that 
the HVI equations ([l])-([3| can be justified, with the ad- 
dition of a nontrivial fluctuation force on the right-hand 
side of ([l]). However, when f2 ^ 1 we are in a quantum 
regime which has seen little experimental exploration, 
and where the vortex equations of motion are rather dif- 
ferent. 

Two key features of the analysis [14] are as follows: 
(i) The widespread assumption of a vortex-quasiparticle 
coupling which is linear in the quasiparticle variables is 
not correct. The vortex is a solitonic excitation of the 
same field as the quasiparticle excitations. Linear cou- 
plings are then forbidden: for the vortex to be a bona 
fide minimum action solution to the equations of motion, 
the lowest-order coupling has to be at least quadratic 
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in all fluctuation variables [11]. One then needs to find 
a "renormalized" coupling to new fluctuation variables 
that are correctly orthogonalized to each other and to 
the vortex "zero mode" ; this turns out to be very compli- 
cated |14| . The renormalized coupling, which is singular 
at low momentum transfer, is indeed quadratic in these 
new variables, (ii) Integration over quasiparticle coor- 
dinates then produces time-retarded, long-range interac- 
tions between different points on a vortex "worldline." 
A nonperturbative path integral treatment (required to 
deal with the singular vortex-quasiparticle interaction) 
then yields extra "memory forces" in the equations of 
motion, which become important in the quantum regime 

n>i. 

(i) Results: At low temperatures, Bose liquids are de- 
scribed by an effective Hamiltonian of the form [TB] : 

H = j d\ (^(^V<i>)2 + Vt?]^ (4) 

with density p — ps + density fluctuations ?7(r), an 
energy functional e[rj, Vrj] whose form depends on which 
superfluid we study, and a superfluid phase <I>(r). This 
Hamiltonian is restricted to lengthscales ^ ao — h/nioCo, 
to energies ^ nioC^, and to velocities ^ Cq, the sound 
velocity, where = Ps{(p£/dr]'^)\,j=o- 

We emphasize the limitations of this hydrodynamic 
formulation. It is valid for both low-density Bose gases 
and dense superfluids like "^He, provided fc^T <C nioC^ 
(so that pn ^ p; in, e.g., ''He superfluid, this means 
T < 0.7 K), and likewise for perturbations of frequency 
ri <C moC^/h. With these restrictions it is valid for arbi- 
trary ratios of the "crossover parameter" = hfl/ksT. 
However it does not include interquasiparticle interac- 
tions, which are very small in this low-T regime and have 
no bearing on the form of the quasiparticle-vortex inter- 
action, but which are essential for the macroscopic trans- 
port of energy and momentum |17j . Nor do we study here 
the role of the boundaries - such geometric effects are 
crucial in understanding the vortex mass [7], although 
they hardly affect the vortex-quasiparticle interaction. 

(a) Equations of Motion: The results are more trans- 
parent when Fourier transformed. Defining R„(f) = 
(ir2Ri,(ri)e*^*, we write the equation of motion as 

RU^) = A,{n,n^)Fj{n), (5) 

where riq is the quasiparticle distribution over momen- 
tum q, and the total "driving force" 

F = Facin) - q,K X 3{n) + F^'"=(r!, nq) (6) 

sums an external driving field, an internal local trans- 
verse force f± = —QuK X J(r2), where J — psVs + PnVn 
is the total current, and a fluctuation term ¥■^^'^'^{^1, Uq). 
The 2x2 "admittance matrix" Aij = A^^Sij where 

a" = D-'^[-n'^M^{n,nq) + inD\\{n,nq)] 

A-^ = D-'^layKpn- aa,\n\d±{n,nq)] (7) 
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FIG. 1: The longitudinal damping coefficient _D||(f2,r). 
Main figure: The dependence of (Cl, T) on Cl, normalized to 
its zero frequency value Do{T). Inset: The coefiicient Do{T), 
proportional to T'*, plotted as a function of the dimensionless 
temperature {To/h)kBT. 

where a are the usual Pauli matrices, and 

D(f7,nq) [n^My - inD\\f - [npQ - i\n\dA_f (8) 

If we write Aij = flaij, then a{fl,T) is analogous to a 
conductivity tensor, with -D||(^l) playing the role of the 
longitudinal resistivity. 

The key difference between ([5|-([8]) and previous re- 
sults for this problem lies in the frequency dependence of 
£'|[(f2,nq), dj_{fl,nq), My(i},nq), and in the correlator 
X»j(r!,nq) = (i^^'"^(l],nq)F/'""(-r!,nq)) of the fluctu- 
ation force F'^'"'^. Quite generally we flnd that when the 
quasiparticles are in equilibrium at temperature T, the 
"viscous" terms [i.e., D||(17,T), dj_{n,T), and Xiji^,T)] 
in (U-Q take the form Qin,T) = /(T)g(O). The effec- 
tive mass AIy{n,T) has a more complicated behavior. 

(b) Quantum- Classical crossover. Consider first 
D\\{Cl,T), shown in Fig. [l] Over the whole range of 
n and T, D||(f2,T) = Do{T)gD{Cl), where Do{T) is 
just the coefficient in ([2|, and where goi^) decreases 
smoothly from (?d(0) = 1 in the classical limit to goi^ ~^ 
oo) = 1/16 in the quantum limit [18[. The fluctuation 
correlator Xiji^jT), shown in Fig. ^ is a little more 
complicated. Ordinarily, one expects the noise correla- 
tor in a quantum Langevin equation to have the form 
Xql{^,T) ^ f{T)tlcoth^fl, i.e., a strictly increasing 
function of Cl. However Fig. shows a quite differ- 

ent behavior: like D^^{(l,T), the longitudinal correlator 
Xij(f2,T) decreases smoothly with il, now with the lim- 
iting behavior 
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known hydrodynamic result |19j 



FIG. 2: The correlator Xij of the fluctuating force, plotted as 
a function of Cl (top) and the dimensionless time kBT{t — s) /h 
(bottom). The plots are normalized to xa a-t zero frequency 
(top) and at zero time (bottom). The arrow at t — s = in the 
bottom figure represents a 5-function contribution (see text). 



where Xo{T) = x^^{^ ^ 0,T) = ADo{T), with A = 
kBT/LzTT (a relationship coming from the fluctuation- 
dissipation theorem). On the other hand, in both limits 
the transverse part £ijkXij of Xiji^^T) is zero; it rises 
to a maximum value ~ 0.2x1 (T) in the crossover regime 
f2 - 1. 

A Fourier transform back to the time domain [see Fig. 
[2|b)] reveals an initial S function in Xii{t~ T) [because 
Xiii^^T) is everywhere finite and positive] followed by 
a slow decay oc (t — s)~^. Similar behavior is found for 
Z)|| (i— s, T) [M] . The transverse correlator rises from zero 
at zero time, and decays more rapidly [like It — s)~'^] at 
long times. 

The transverse function d^^{yi,T) also shows a charac- 
teristic crossover behavior - we do not elaborate here 
because d±{Vt,T) is always very small compared to 
D\\(yi^T). Full details on all these functions are found 
in the Supplemental Material [14] . 

Finally, consider the vortex inertial mass Mu(r2, T) ap- 
pearing in ([5|-([8]). This is well-known to depend on the 
sample geometry [7]; for a circular container of radius 
Ro S> cto, in the J7 — > limit, we easily verify the well- 
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+ 7is + 1/4 



(10) 



where 7^ is Euler's constant. Naively, one expects that 
in the quantum limit 3> 1, the ln(i?o/ao) factor in ( 10 ) 
will be replaced by hi[co/ao^) (the length scale Co/flo be- 
ing set by the distance quasiparticles can travel in a time 
~ il^^). However the actual behavior is more subtle: 
there is a "radiation reaction" term cx d?'R,y/dt^ in the 
equation of motion, analogous to that in electrodynamics, 
and to deal with this one must go beyond the expansion 
in powers of Ri,/co being used here. This problem lies 
outside the scope of the present paper [20| . 

(c) Real-time dynamics: Remarkably, the results given 
above allow us to write simplified local equations of mo- 
tion in both quantum and classical limits. In the classical 
limit, Fourier transforming back gives precisely the HVI 
equations ([l])-([3]), but with an added noise term: 

M,R„ - /m - fqp - ^ac[t) = F^fi(t) (11) 

where the classical noise force Fyj'j^(t) has the correlator 



xg')(t-5,T)^x!l(r)M(^-.) 



(12) 



i.e. an entirely longitudinal Markovian noise. However 
this equation is only meaningful on coarse-grained time 
scales ^ h/ksT] for shorter times, the time- retarded 
nature of the correlations becomes crucial, and as Fig. [2] 
makes clear, the fluctuation correlator Xij 

{t - s, T) then 

becomes anisotropic and highly non-Markovian, and the 
HVI equations simply do not apply. 

In the opposite quantum limit 3> 1, one may 
again write a local equation like (11), of HVI form, 
again with an added noisy fluctuating force. However, 
now the coefficients are different; Do{T) in ^ IS re- 
placed by Z?o(T)/16, and the quantum noise correlator 
X^\t -s,T) = ^xif (t - s,T) (again entirely Ion- 
gitudinal). In this limit, valid for time scales ^ h/ksT 
(but » h/moCg), the coefficients in these "quantum HVI" 
equations arise solely from the (5-function contributions 
to the correlation functions - all retarded parts are sup- 
pressed. 

We would like to emphasize how unusual these results 
are in detail. It is quite remarkable to have 2 equations 
of exactly the same form (but quite different coefficients) 
in these two limits, but with a quite different form in 
the crossover between them; it is more illuminating to 
look at it in frequency space, as above. And yet, very 
surprisingly (at least to us), the lordanskii force is quite 
unaffected by this - apart from the very small correction 
term d±{n, T), the lordanskii force is independent of fre- 
quency, and can be treated as entirely local, and ([3| is 
reproduced exactly in our derivation (which is quite dif- 
ferent from previous scattering theory calculations). 
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(ii) Experimental Implications: The results above 
justify the HVI equations ^SjiOj, and the phenomenology 
based on these [10] , in the classical regime. However away 
from this regime we find a quite new phenomenology. 
Clean experimental tests of this will require (a) that the 
vortex not be coupled to some other object (e.g., charged 
ions), which change its natural dynamics, and (b) that 
the vortex be coupled to the natural excitations of the 
system (as opposed to, e.g., a source of external quasi- 
particles, which can couple linearly to the vortex). The 
results will also change in situations where the vortex is 
being "dragged" by some external time- varying potential 
[5], since such potentials may strongly distort the vortex 
in the region where they act. 

The most obvious direct experimental realization of the 
results here would be in 2-dimensional Bose-condensed 
atomic gases; the dynamics of single vortices can then be 
tracked in experiments [21j (e.g., in their spiraling out 
from the trap center), and the viscous coefficients can 
then in principle be extracted from such measurements 
[52]. A detailed treatment using the present equations 
is quite lengthy, and will be presented elsewhere. Our 
results are also clearly relevant to experiments on tur- 
bulence in superfluid "^He, which have recently begun to 
probe the quantum regime [23 , and any theory of vortex 
tunneling in ^He or cold gases must include the viscous 
effects described here (which are very far from being de- 
scribed by a simple Caldeira-Leggett coupling |24J). One 
difficulty with turbulence or tunneling is that in most ex- 
perimental cases the vortices are 3-dimensional objects, 
and the vortex line may distort in many ways that are 
impossible to capture analytically. It would nevertheless 
be interesting to extend, at least numerically, the existing 
theories of quantum turbulence [25' and vortex tunneling 
[26] in ^He to include the effects discussed here. Finally, 
we emphasize that the results given here are not applica- 
ble to fermionic superfluids like '^He or superconductors 
- the form of the vortex-quasiparticle interaction is quite 
different in these systems. 

To summarize: within the constraints of the low-T hy- 
drodynamic picture, we find that the HVI equations can 
be justified in a purely quantum-mechanical treatment, 
with the addition of a fluctuation noise term, provided 
one is in the classical regime (l <^ 1. Outside this regime 
one needs to use a more general set of equations, which 
show strong memory effects in the time domain. 
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